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Abstract 

The  deformation  in  a  rate-dependent  solid  containing  a  pair  of  closely- spaced  spherical 
voids  is  analytically  determined.  Results  are  obtained  which  are  valid  in  the  material  that 
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September,  1988 


Introduction 

In  this  paper,  the  finite  deformation  and  stress  fields  in  the  vicinity  of  a  pair  of  spherical 
cavities  in  an  all-around  infinite  rate-dependent  solid  are  examined.  The  cavities  remain 
traction-free  as  a  uniform  stress  state  is  remotely  applied  to  the  surrounding  material.  The 
separation  of  the  cavities  is  assumed  to  be  much  smaller  than  their  radii.  This  assumption 
enables  the  use  of  a  theory  of  plane  stress  to  describe  the  deformation  of  the  thin  piece  of 
material  which  separates  the  cavities.  The  theory  used  here  is  similar  to  that  described  by 
Adkins,  et  al  (1954)  for  hyperelastic  solids. 

The  problem  dealt  with  here  is  of  potential  importance  in  connection  with  the  rupture 
and  failure  of  solids  by  the  micromechanisms  of  void  growth  and  coalescence.  Thus,  an 
understanding  of  the  interaction  of  a  pair  of  closely-spaced  cavities  serves  to  shed  light  on 
a  possible  crack-advance  process.  Such  results  might  also  eventually  serve  as  useful  guide 
to  the  development  of  constitutive  models  that  accurately  represent  the  behavior  of  porous 
solids  at  high  levels  of  porosity.  Existing  models,  such  as  that  proposed  by  Gurson  (1977) 
were  designed  to  model  materials  that  contain  a  sparse  population  of  voids. 

Within  the  linearized  theory  of  elasticity,  the  problem  of  the  interaction  between  a 
pair  of  spherical  cavities  in  an  isotropic,  homogeneous  solid  was  solved  by  Sternberg  and 
Sadowsky  (1952).  They  obtained  solutions  to  the  global  problem  in  the  form  of  infinite 
series,  and  did  not  need  to  assume  that  the  voids  were  closely-spaced. 

In  the  context  of  finite  elasticity,  single-void  problems  have  been  examined  by  several 
authors,  including  Green  and  Shield  (1950),  Ball  (19S2),  and  Horgan  and  Abcyaratne 
(1986).  Single  cavity  problems  in  rate-dependent  solids  were  examined  by  Abeyaratne 
(1988),  who  was  particulary  interested  in  the  prediction  of  void  nucleation.  An  analysis 
pertaining  to  contiguous  voids  in  a  class  of  hyperelastic  solids  was  also  recently  done  (Blume 
(1988)). 

In  the  first  section  of  this  paper,  some  preliminaries  on  finite  deformations  of  rate 
dependent  solids  arc  assembled.  The  particular  constitutive  class  which  will  be  used  for 
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the  boundary-value  problem  at  hand  is  introduced  and  discussed.  In  Section  2,  the  analysis 
of  the  two-cavity  problem  is  carried  out,  and  in  the  last  section,  the  results  are  discussed. 


1.  Notation,  Preliminaries  on  Finite  Deformations  of  Rate-Dependent  Solids 
Throughout  this  paper,  E3  stands  for  three-dimensional  Euclidean  point  space.  Lower  case 
and  capital  letters  in  boldface  designate  vectors  and  (second-order)  tensors,  respectively. 
The  same  letter  in  lightface  appearing  with  one  or  more  subscripts  will  signify  the  appro¬ 
priate  components  of  the  vector  or  tensor  in  either  cylindrical  or  Cartesian  coordinates. 
If  A  is  region  in  E3 ,  /  £  Cm(A)  (m=l,2,3,...)  signifies  that  /  is  m  times  continuously 
differentiable  on  A. 

If  1Z  is  the  region  in  E3  occupied  by  a  solid  in  its  reference  configuration,  a  motion  is 
a  mapping  y  described  by 

y  =  y(x,f)  =  x  +  u(x,f)  on  TZ,  {t  >  0),  (1.1) 

where  x  is  the  position  vector  of  a  generic  point  in  7Z,  y(x,  t )  is  the  position  vector  of  its 
deformation  image  at  time  t,  u  is  the  associated  displacement  field,  and  t  is  time.  The 
region  occupied  by  the  body  at  time  t  will  be  called  7£(,  and  at  each  t,  the  mapping  from 
TZ  to  lZt  is  assumed  to  be  one-to-one.  The  motion  y  of  7Z  is  taken  to  be  twice  jointly 
continuously  differentiable  with  respect  to  position  and  time. 

The  solids  which  will  be  dealt  with  here  are  incompressible  and  can  thus  sustain  only 
locally  volume-preserving  deformations.  The  deformation  (1.1)  is  locally- volume  preserving 
if  and  only 

J  =  detF  =  1  on  7?.,  ( t  >0),  F  =  Vy.  (1-2) 

Here,  F  is  the  deformation  gradient  tensor  and  J  is  the  Jacobian.  The  letters  C  and  G 
stand  for  the  right  and  left  Cauchy-Green  strain-tensor  fields: 

C  =  FrF,  G  =  FFronft.  (1.3)' 

Both  C  and  G  are  symmetric,  positive-definite  tensor  fields  with  common  (positive)  prin¬ 
cipal  values.  If  D  denotes  the  rate-of- deformation  tensor,  then 

D(y,f)  =  symVv(y,f)  on  lZt,  t>  0  (1.4) 

1  The  superscript  T  indicates  transposition. 
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in  which  v  is  the  Eulerian  velocity  field: 

v(y,t)  =y(x(y,t),t)  on7Zt,  t>  0.  (1-5) 

Here  the  dot  denotes  the  derivative  with  respect  to  time  (x  fixed),  and  x(y,t)  is  the  inverse 
of  y(x,  t)  at  time  t:  x(y(x,<),<)  =  x  for  all  x  in  71.  The  rate  of  deformation  tensor  is  thus 
expressible  as 

D(y(x,t),<)=^{F(x,i)F-‘(x,()  +  F-T(x,()FT(x,()}.  (1.6) 

If  a  is  the  Piola  stress  tensor  accompanying  the  deformation  at  hand,  equilibrium,  in 
absence  of  body  forces,  demands  that 

diva  =  0,  aFr  =  F <rT  on  71.  (1.7) 

In  terms  of  the  Cauchy  stress  field  r,  which  is  linked  to  the  Piola  stress  tensor  through 

r( y,  t)  =  a(x,  t)Fr(x, <),  y  =  y(x),  (1.8) 

equilibrium  may  be  written  as 

divr(y,  t)  =  0,  r(y,  t)  =  rT(y,  t)  for  all  y  6  7lt  (1.9) 

Note  that  inertial  effects  are  neglected,  in  the  analysis  carried  out  here. 

The  stress  response  of  materials  at  hand  is  assumed  to  be  well  described  by  a  purely 
rate-dependent,  power-law  type  constitutive  law.  In  particular,  if  u  and  7  are  material 
constants  with  v  >  0,  then 

r  =  —pi  -f  1  /(PD  on  7vt.  (1.10) 

where  p  =  p(y,  t)  is  the  pressure  field  which  accommodates  the  constraint  of  incompress¬ 
ibility,  1  is  the  identity  tensor,  and 


<1  =  | D 1  =  \J D ij D ij  on  Tit. 


(1.11) 


Note  that  the  second  of  (1.9)  is  satisfied  automatically.  On  account  of  (1.8),  (1.3),  and 
(1.6),  the  Piola  stress  for  such  materials  may  be  written  as 


a  =  -PF~T  +  Kd?  {fC"1  -  (F-T)  }  on  K. 


(1.12) 


In  order  to  understand  the  behavior  of  such  a  solid,  consider  a  state  of  uniaxial  stress. 
If  the  body  is  in  a  state  of  uniaxial  stress,  then  only  one  stress  component,  say  rn  is 
nonzero  and  depends  only  on  time;  all  others  vanish  identically  in  7£.  One  motion  and 
pressure  field  which  corresponds  to  such  a  state  is 


yi(x,t)  =  A(t)x1,  y2(x,t)  =  \  1/2(f)x2,  y3(x,/)^A  1/2{t)x3 


1/2 /a  V+' 


(1.13) 


Here  A  >  0  for  t  >  0,  and  A  is  a  twice  continuously  differentiable  function  of  time.  One 


rn  =  v 


'-(!) 


7/2  +  1 


Further, 


U‘  /rn(s)\^  ) 

^ — i — J  ds>,  c  constant. 


For  Tn(t)  =  Tt,  (r  constant) 


A(t)=cexpjGr+1 


(1.14) 


(1.15) 


(1.16) 


The  response  is  shown  in  Figure  1. 


2.  The  Two-Cavity  Problem 


If  {ei,e2,e3}  are  base  vectors  of  a  Cartesian  coordinate  frame,  the  solid  initially  occupies 
the  unbounded  region  7?.,  described  by 

H  =  {x|  |x  -  <5e3 j  >1,  |x  4-  6e3|  >  1}.  (2.1) 

The  coordinate  system  is  chosen  so  that  the  holes  have  unit  radii;  their  centers  are  located 
at  x3  =  ±6,  with  6  >  1.  The  voids  are  thus  separated  by  a  distance  of  2(6  —  1).  In  what 
follows,  it  is  assumed  that 

1  <  6  <  2;  (2.2) 

this  assures  that  the  cavities  are  closely- spaced.  In  order  that  the  cavity  surfaces  remain 
traction  free,  the  nominal  traction  vector  must  vanish  there: 

on  =  0  on  671  ~  {x|  |x  ±  6e3|  =  1},  (2.3) 

or  equivalently,  the  Cauchy  traction  vanishes  on  the  deformed  cavity  walls,  whence 

Tn  =  0  on  0Rt.  (2.4) 

In  the  preceding  expressions,  n  is  the  outer  unit  normal  vector  to  dl v;  n  is  the  outer  unit 
normal  vector  to  the  cavity  surfaces  in  the  current  configuration.  At  large  distances  from 
the  voids,  the  solid  is  subjected  to  a  uniform  state  of  stress  which  has  axial  symmetry 
about  the  x3-axis  as  well  as  symmetry  about  the  plane  .r3  =  0.  The  problem  geometry  is 
shown  in  Figure  2. 

The  main  purpose  of  this  analysis  is  to  characterize  the  deformation  in  the  material 
which  separates  the  voids.  For  moderate  values  of  the  remotely  applied  stress  field,  and  for 
large  times,  it  appears  reasonable  to  assume  that  the  largest  strains  and  strain-rates  occur 
in  this  region,  and  that  rest  of  the  solid  might  well  undergo  infinitesimal  deformations  at 
low  rates. 


Thus,  consider  the  region  71,,  that  separates  the  voids  in  the  undeformed  configuration 
(Figure  3): 

71,  =  {x|  r  <  a,  —  h(r)  <  x3  <  h(r)} ,  r  =  yj x\  +  x\,  h(r )  =  S  -  \/l  —  r2,  (2.5) 

and  2(6  —  l)<a<lisa  real  number.  Assume  also  that  the  (axially  symmetric)  remote 
load  in  the  global  two-cavity  problem  produces  normal  and  shear  traction  distributions 
crn(x3,t)  and  ca(x3 ,t)  on  the  lateral  boundary  of  the  disk  71,,  so  that  the  components  of 
the  Piola  stress  in  cylindrical  coordinates  (r,d,x3)  obey 


0rr(r,x3,t)  =  an(x3,t), 
* 

°r3  (r,x3,t)  =  at(x3,t), 


r  =  a,  —h(a)  <  x3  <  h(a),  t  >  0. 


The  symmetry  of  the  problem  mandates  that 

Zn(x3,i)  =  crn(-x3,t)  as(x3,t)  =  -as(-x3,t),  -h(a)  <  x3  <  h(a).  (2.7) 

The  faces  of  the  disk  i3  =  ±h(r)  (0  <  r  <  a)  are  traction-free.  In  terms  of  the 
cylindrical  components  of  Piola  stress,  this  last  condition  is  equivalent  to 


±  arr(r,x3,t)h'(r)  -  (jr3(r,x3,t)  =  0, 
±  <*3r (r,x3,t)h'(r)  -  (j33(r,x3,t)  =  0 


on  x3  =  ±h(r),  h'  =  — ,  (0  <  r  <  a).  (2.S) 


At  this  point,  it  is  useful  to  recall  that  for  axially  symmetric  motions,  the  equilibrium 
equations  (1.7)  in  terms  of  the  cylindrical  coordinates  ( r,0,x3 )  reduce  to 

d  l  .  d 

o  ^rr  4"  (^rr  GOB)  d"  „  O r3  ~~  0, 

dr  r  9X3  (2  9) 

d  d  i  1  j 

■X“CT3r  d-  -X - <7.33  d - <7r3  —  U. 

or  or  3  r 

Using  plane-stress  approximative  assumptions,  the  motion  in  the  disk  7v*  induced  by 
the  edge-load  distributions  a„  and  a,  will  bo  calculated.  Owing  to  the  symmetries  inherent 
in  the  problem,  any  motion  y  of  7 Z  (and  hence  also  of  the  subregion  7v»  of  7v)  has  the  form 


ya(x,t)  =  R(r,x3,t)xn/r,  y3(r,x3,t.)  =  Z(r,.r3,t), 
R(r,x3,t)  =  R(r,  —x3,t),  Z(r,  x3,t)  =  -Z(r,-x3,i), 


(a  =  1,2)  (2.10) 
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with  R  =  0  when  r  =  0,  but  R  and  Z  nonvanishing  and  twice  continuously  differentiable 
for  (0  <  r  <  a).  Further,  the  nominal  stress  components  obey 


orr(r,x3,t)  =  orr(r,-x3,t),  agg(r,x3,t )  =  aee(r,  - x3 ,  t),  a33(r,x3,t)  =  a33(r,  -x3,  t), 
Or 3(r,x3,t)  =  -ar3 (r,-x3,t),  a3r(r,x3,t)  =  -a3r(r,-x3,t)i 

(2.11) 

all  other  stress  components  vanish. 

Let  II  stand  for  the  midplane  (13  =  0)  of  7£„,  and  if  /  is  a  function  defined  on  7 £»,  / 
will  be  used  to  denote  its  restriction  to  II.  The  aim  of  the  following  is  to  characterize  the 

o  o 

deformation  parameters  R  and  dZ(r,0,t)/dx3  =  A (r,  t),  along  with  the  stress  components 
ffrr  and  age  by  assuming  that 


cr33  =  0  on  II, 


—  agg  —  0  on  71,. 
(  3 


(2.12) 


These  may  be  justified  on  the  basis  of  the  geometry  of  the  problem,  the  second  of  (2.8), 
and  the  second  of  (2.9). 

Integrating  the  first  of  the  equilibrium  equations  (2.9)  across  the  thickness  of  TZ ,  leads 


to 


r{r)  fo  1  1 

/  <  ~^°rr  +  -  (arr  -  agg)  ^  dx3  +  ar3 

J-h(r)  {or  r  J 

The  first  of  (2.8)  and  the  assumptions  (2.12)  enable 


r3  =  /i(r) 


=  0. 


(2.13) 


z3=-h(r) 


~  Ogg)  + 


=  0  on  n. 


(2.14) 


The  boundary  condition  in  the  first  of  (2.C)  corresponds  under  the  present  circumstance 
to  the  stipulation  that 

arr(a,t)  =  o„(a,Q,l).  (2.15) 

Note  that  the  integral  of  the  shear  traction  as  across  the  thickness  vanishes  by  virtue  of 
the  symmetry  condition  that  appears  in  the  second  of  (2.7). 

With  a  view  towards  expressing  the  stresses  arr  and  dgg  in  terms  of  the  deformation 

©  o 

parameters  R  and  A,  note  that  (2.10)  implies  that  the  nonzero  cylindrical  components  of 


the  deformation  gradient  tensor  F  are 

Frr  =  R,r,  F$e=R/r ,  F33  =  A  =  8Z(r,  0,  t)/dxi  on  II. 
Here  and  in  what  follows,  dR/dr  =  ilr,  while  dR/dt  =  R<t. 

o 

Incompressibility  demands  that  det  F  =  1  on  II,  whence 


(2.16) 


A  = 


on  n. 


R.rR 


The  stress  components  are  now  found  from  (1.12)  to  be 


+  1  J-,  aee  =  \  -p  +  udy^-  1  4, 

R.r  I  R.r  I 


RR,r 


i  —p — vd1 

o  o 

R.t  R,rt 

O  "t“  O 

l 

R  R,r 

and 


d  = 


R,»  ( kX  R,rtk.,  . 

2  1  -5 —  +  ~  H - o — —  on  II. 

R,r  \R  )  R,rR 


(2.17) 


(2.18) 


(2.19) 


The  last  of  (2.18),  in  conjunction  with  the  plane  stress  assumption  appearing  in  the  first 
of  (2.12)  gives 


p  =  —  vd? 


R.t  ,  R,rt 

®  i  ° 

R  R,r 


on  n, 


whence 


O  rr  =  V(P 


J 

+ 

to 

j3j» 

_ 1 

dee  =  vd1 

0R,t  R.n 

—  o  1  o 

1 

W 

fcf 

a* 

R.r 

R  R,r 

2_ 

R' 


(2.20) 


(2.21) 


On  entering  (2.20)  into  (2.14),  one  is  led  to  a  nonlinear,  homogeneous  partial  differential 

o 

equation  for  R  with  nonconstant  coefficients.  This  equation  involves  the  second  derivative 

o 

of  R  with  respect  to  r,  and  its  first  derivative  with  respect  to  time.  The  boundary  condition 

o 

(2.15)  requires  that  R  conform  to 


udy 

Xr 


R  "R.r 


=  0n(a,O,t). 


(2.22) 
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Of  particular  physical  interest  in  the  problem  at  hand  is  the  large-time  behavior  of 

o 

R.  This  is  perhaps  when  the  neglection  of  inertial  effects  is  most  valid.  It  is  reasonable  to 

o 

assume  that  R  is  of  the  form 


R(r,t )  =  A(t)q(r),  0  <  r  <  a,  t  >  0, 
with  A  and  q  twice  differentiable.  In  this  case,  (2.19)  gives 


and  (2.21)  leads  to 


d  =  V6— ,  A 
A 


dA 
dt  ' 


(2.23) 


(2.24) 


in  which 


°rr (r,t) 


m(t) 

q'(r)' 


„  m(t)  _  , 

1 796  =  r——  on  n,  q 
q(r) 


dq 
dr ' 


m  =  2uQ>’1^2 


(2.25) 


(2.26) 


Substitution  from  (2.25)  into  the  equilibrium  equation  (2.14)  and  rearranging  leads  to  the 
following  ordinary  differential  equation  for  q: 


If  one  sets 


~r  [<l<l"  +  (<?')2]  +  W  +  r^qq'  =  0, 


0  <  r  <  a. 


(2.27) 


v(r)  =  q{r)q'(r),  0  <  r  <  a, 


then  (2.27)  reduces  to  the  first  order  equation 


v 


i 


(2.2S) 


(2.29) 


the  solution  to  which  is  given  by 


v(r )  =  crh(r). 


(2.30) 


Here,  c  is  a  constant.  Bearing  in  mind  that  R(0,t)  -=  0,  one  finds  on  integrating  (2.2S) 


that 


(2.31) 
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o 

Appealing  to  the  definition  of  h  in  the  last  of  (2.5)  and  the  expression  for  R  in  (2.23),  one 


o 

Plots  of  R  versus  r  for  fixed  t  are  shown  in  Figure  4. 

The  function  A(t)  is  related  to  an(a,0,t)  through  (2.22),  which  now  gives 


(2.32) 


A(t)  =  expj  J  <ryn+  (a,o,s)ds| 


q'(a )  =  2a- 


6  —  Vl  —  a 2 


(2.33) 


\J8a2  +  1(1  -  «2)3/2  ~  f 

If  an  is  time- independent,  A(t)  is  an  exponential  function  of  time,  while  if  an  is  a  linear 
function  of  time,  A(t)  is  a  constant  multiple  of  the  function  A(<)  in  (1.16),  which  is  plotted 
in  Figure  1. 


According  to  (2.32),  (2.17),  the  out-of- plane  stretch  A  is  given  by 


A(r,  <)  =  A  2(t ) 


6  —  v/l  —  r2 


(2.34) 


which,  on  account  of  the  last  of  (2.5)  may  be  written  as 

kr’V  =  #mry  (235) 

The  boundary  13  =  h(r)  is  thus  approximately  mapped  to  the  plane  y3  =  Z  =  A~2{t)  and 
the  material  in  between  the  spheres  is  flattening.  Figure  5  shows  A  versus  r,  while  Figure 
6  shows  the  deformed  disks  for  various  values  of  A(t)  at  fixed  t.  In  particular,  note  that 


A(0,<)  = 


(2.36) 


v~’  '  A2(t)(6-lY 

o 

is  the  maximum  value  of  A  at  fixed  f;  the  minimum  occurs  on  the  edge  of  the  disk  (r  =  a). 
As  far  as  the  stresses  are  concerned,  the  Nominal  stresses  and  Cauchy  stresses  are 

found  from  (2.32),  (2.21),  (1.8): 

.  _  ,  .  Sr  —  r\J  1  —  r2  „  m(t)r 


brr  —  m(t)- 


vee  = 


^r2  +  |(1_r2)3/2_|  ^r2  +  |(1_r2)3/2_|  (2.37) 


Trr  =Tge  =  ™(t), 

with  m  given  by  (2.26),  (2.32).  Note  that  the  Cauchy  stresses  are  constant. 
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3.  Discussion 


In  this  analysis,  it  was  assumed  that  the  material  between  the  voids  does  not  rupture.  This 
was  reflected  in  the  stipulation  that  R  vanishes  at  r  =  0.  If  this  requirement  is  removed 
and  replaced  with  R  >  0  at  r  =  0,  the  possibility  of  a  hole  opening  at  the  center  of  the 
disk  is  admitted.  In  this  case,  in  order  that  the  newly-formed  hole  be  traction-free,  one 
demands  that 

oTT  =  <73r  =  0  if  r  =  0.  (3-1) 


Pursuing  the  plane  stress  analysis  exactly  as  before  leads  one  to 


R(r,t)  =  A(t)^6,-*  +  j  (1  -  r*)3/2  - 


(3.2) 


with  A  >  0  a  function  of  time  and  b  constant,  b  <  2/3.  The  boundary  condition  (2. 22^ 

o 

must  hold  again  in  this  case,  but  in  place  of  the  requirement  that  I?(0,  t)  =  0,  one  has  from 
the  first  of  (3.1)  that 


<7rr(0,t)  = 


vdi 

o 

R,r 


R,t 

R  R,r 


0. 


(3.3) 


r=0 


The  second  of  (3.1)  is  satisfied  in  average  across  the  thickness  by  virtue  of  symmetry. 
However,  the  condition  (3.3)  can  be  met  only  if  b  =  0  so  that  plane  stress  deformations  of 
the  form  (2.23)  with  holes  emerging  at  the  origin  are  not  present. 

Although  no  restrictions  (other  than  symmetry  conditions)  were  put  on  the  remote 
stress  field  in  the  global  problem,  the  plane  stress  assumptions  are  perhaps  more  believable 
in  situations  in  which  the  disk  71*  is  itself  in  a  tensile  field  (so  that  it  is  thinning ); 
presumably  this  situation  would  correspond  to  remote  tensile  loading. 

The  validity  of  the  plane  stress  assumptions  is  increased  if  boundary  of  the  disk  7v, 
is  chosen  so  that  the  magnitude  of  the  shear  traction  aa(x3,t)  is  minimized  across  the 
thickness  of  71 *.  As  yet,  a  means  of  connecting  the  tractions  a„  and  aa  with  the  remote 
loading  is  unavailable;  this  is  the  subject  of  future  work.  The  effects  of  inertia  on  the 
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deformation  of  7£,  Eire  also  being  examined;  plane  stress  may  not  be  a  valid  assumption  in 
the  dynamic  case. 
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Figure  Captions 


Figure  1:  Response  of  the  power  law  viscous  solids  in  uniaxial  stress; 

A.(i)  =  A(()exp{-(I),/,''+,,l±l}. 

Figure  2:  Problem  geometry. 

Figure  3:  The  region  TZt . 

o 

Figure  4:  R  versus  r. 

Figure  5:  The  out-of-plane  stretch  A  versus  r. 

Figure  6:  Deformed  disks. 
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